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METHOD AND APPARATUS FOR TWO DIMENSIONAL SURFACE PROPERTY 
ANALYSIS BASED ON BOUNDARY MEASUREMENT 

CONTRACTUAL ORIGIN OF THE INVENTION 
[0001] This invention was made with United States Government support under Contract 
No. DE-AC07-99ID 13727 awarded by the United States Department of Energy. The United States 
Government has certain rights in the invention. 

BACKGROUND OF THE INVENTION 

[0002] Field of the Invention: The present invention relates to nondestructive material 
evaluation and, in particular, to a method for characterizing and evaluating integrity and physical 
properties of electrically conducting material structures. 

[0003] Description of Related Art: Structures and the materials comprising those 
structures are often evaluated for integrity and other physical properties using a variety of 
nondestructive evaluation techniques. These techniques include thermographic, optical, acoustic, 
radiographic (e.g., x-ray), and electromagnetic procedures. 

[0004] Electrical impedance measurement and analysis of structures and their surfaces is 
of particular interest because often, the impedance characteristics of a surface may correlate to other 
physical phenomena, such as strain on the surface, thermal characteristics of the surface, 
photosensitivity of the surface, and physical integrity of the surface. 

[0005] In the area of strain measurements, conventional foil type electrical strain gauges 
are often used to measure deformation of a structure or deformation of the surface of a structure. 
However, strain gauges typically cover only a small area of the structure. As a result, strain gauges 
are typically used in specific areas of a structure where a problem, such as stress concentration or 
orientation of concern, is likely to exist or where measurements are desired due to the geometry and 
configuration of the structure at that specific area. The cost and effort required to cover a significant 
surface area of a structure with resistance strain gauges may be prohibitive. Additionally, a 
technique using multiple strain gauges would likely require the presence of conductive, such as 
wire, leads extending to each strain gauge distributed across the surface to be observed, which may 
be impractical in many applications. 
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[0006] Techniques exist for measuring and modeling physical properties of the surface 
of a structure, when the surface is electrically conductive. See, for example, U. S. Pat. No. 
5,165,794 to Ortiz for a Method For The Thermal Characterization, Visualization, and Integrity 
Evaluation of Conducting Material Samples or Complex Structures. However, like a strain gauge 
implementation, the Ortiz patent requires wire leads attached to measurement points distributed 
throughout the surface of the structure, which may be impractical in many applications. 

[0007] Many techniques, such as electrical impedance tomography, exist for measuring 
and analyzing electrical impedance characteristics of an object by only placing measurement probes 
around the periphery of the object. Unfortunately, these solutions are typically volumetric solutions 
requiring the object to be somewhat conductive throughout the volume of the object. Additionally, 
tomography techniques are generally used to characterize resistivity through the object as a whole in 
three dimensions, or at least of a cross section through the object. Thus, for analysis concerned with 
the surface characteristics of an object, these electrical impedance tomography techniques are 
generally inadequate. 

[0008] A method and apparatus is needed to nondestructively measure and analyze 
electrical resistance parameters of a conductive film, or a conductive surface of a structure, while 
only requiring the use of measurement points around the periphery of the conductive film or 
conductive surface. Additionally, a method and apparatus is needed to extrapolate the electrical 
resistive parameters of the structure to other physical properties such as thermal characteristics of 
the structure at the surface, photosensitivity of the surface of the structure, strain on the structure 
surface, and physical integrity of the structure surface. 

BRIEF SUMMARY OF THE INVENTION 
[0009] One embodiment of the present invention includes a method for determining 
physical properties of a conductive film. A plurality of probe locations may be selected around a 
periphery of the conductive film. A plurality of measurement lines may be defined as the line 
segments that may be drawn between each probe location and all other probe locations in the 
plurality of probe locations. Electrical resistance may be measured along each of the plurality of 
measurement lines. Based on the measured electrical resistance values, a lumped parameter 
resistance model may be developed with lumped resistance values existing along each measurement 
line. Using linear algebra, the lumped parameter resistance model may be determined based on the 
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measured resistance values. Using the lumped parameter resistance model, electrical resistivity may 
be modeled at any selected location encompassed by the plurality of probe locations. If the 
conductive film includes a correlation between electrical resistivity and another physical property, 
the correlation may be used to extrapolate from the modeled resistivity at the selected location to the 
desired physical property at the same selected location. A profile of the conductive film may be 
developed by determining the resistivity and other physical properties at a plurality of selected 
locations. 

[0010] In another embodiment of the present invention, the conductive film to be 
analyzed may be applied to the surface of a structure. Once the conductive film is applied, the 
conductive film may be analyzed as described to arrive at resistivity or other physical property at the 
selected location. An electrical resistivity value may be estimated at a selected location on the 
conductive film encompassed by the plurality of probe locations. The structure may be an 
essentially two-dimensional object on which the surface is to be analyzed. Additionally, the 
structure may be a three-dimensional structure, wherein the conductive film covers a surface or 
multiple surfaces of the structure. Additionally, the surface(s) covered may be curved. In other 
words, the surface(s) may be nonplanar. 

[0011] In another embodiment of the present invention, the previously described analysis 
may be performed on a structure that already incorporates the conductive film. 

[0012] Yet another embodiment of the present invention includes a system configured 
for determining surface properties of a structure bearing a conductive film over a surface of the 
structure. The system includes a plurality of probes adapted for of measuring an electrical resistance 
when the plurality of probes are placed at the plurality of probe locations around the periphery of 
the conductive film. The plurality of probes connect to a signal controller. The signal controller may 
be configured to select any pair of the plurality of probes at any given time for measuring the 
resistance values along all the measurement lines. A signal sampler, operably coupled to the signal 
controller, may be configured to sample the electrical resistance between the currently selected pair 
of the plurality of probes. The sample may be digitized and sent to a suitably programmed 
processor, which may then be used to perform the method described above to determine resistivity 
or other physical properties at a selected location. By analyzing a plurality of selected locations, the 
system may be used to develop a profile of the resistivity or other physical properties across the 
surface of the conductive film or surface of the structure. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0013] In the drawings, which illustrate what is currently considered to be the best mode 
for carrying out the invention: 

[0014] FIG. 1 is a diagram of a circular conductive film indicating probe locations and 
measurement lines on the conductive film; 

[0015] FIG. 2 is a diagram indicating measurement lines and various distances used 
calculating weighted averages used in determining resistivity at a selected location; 

[0016] FIG. 3 is a diagram depicting various analysis triangles that may be used in a 
weighted average analysis to determine resistivity at a selected location; 

[0017] FIG. 4A is a diagram of a rectangular conductive film indicating probe locations 
and measurement lines on the conductive film; 

[0018] FIG. 4B is a representation of a cylindrical structure for acceptance of a 
conductive film; 

[0019] FIG. 4C is a representation of the cylindrical structure with the conductive film 
applied to the surface of the structure; and 

[0020] FIG. 5 is a block diagram of a system for analyzing surface properties of a 
structure bearing a conductive film. 

DETAILED DESCRIPTION OF THE INVENTION 
[0021] An exemplary embodiment of the present invention includes a method of 
analyzing resistance, or other physical properties, of a conductive film 1 10. As shown in FIG. 1, a 
plurality of probe locations 120 may be selected around the periphery 130 of the conductive film 
1 10. The probe locations 120 may be selected in a manner allowing one or more areas of interest, 
where property analysis is desired to be encompassed within a boundary defined by line segments 
extending between each of the adjacent probe locations 120. For example, in FIG. 1, with N(six) 
probe locations 120, the boundary is defined by six line segments (i.e., line segment P1-P2, line 
segment P2-P3, line segment P3-P4, line segment P4-P5, line segment P5-P6, and line segment P6- 
Pl). 

[0022] As a result of selection of the probe locations 120, an analysis mesh is defined by 
line segments, also referred to as measurement lines 150, between each probe location and all other 
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probe locations 120 in the plurality of probe locations 120. In general, N probe locations 120 define 
a total of N*N measurement lines 150 between the probe points at which discrete resistance 
measurements may be taken. Obviously, this total of N*N measurements includes measurements of 
resistivity between a probe location and itself, which is a point rather than a line segment, and 
therefore need not be measured. Additionally, a resistance measurement along, for example, line 
segment P1-P5 will be the same when measured from PI to P5 and from P5 to PI. So, even though 
it is convenient to discuss, and populate a matrix with, N*N discrete measurement values, only 
(N*(N-l))/2 actual discrete measurements are required. 

[0023] A mesh, in the context of this analysis, refers to the resultant set of lines as 
defined above. It is important to note that this mesh does not necessarily refer to the intersections 
between lines as may be thought of in the case of a typical mesh. The analysis method used in the 
present invention does not rely on determining property values at intersection points of 
measurement lines 150. Rather, it relies on measurements around the periphery 130 combined with 
analysis and estimation using these periphery measurements without requiring determination of 
properties at internal intersection points. 

[0024] Resistance measurements are obtained along each of the measurement lines 150 
to arrive at a set of measured resistance values 140. These measured resistance values 140 are 
shown in the figures with a capital R followed by the points defining the line segment along which 
the measurement is taken. Resistance may be measured using various techniques well knowri in the 
art. For example, using a simple Ohmmeter to directly measure resistance, placing a specific voltage 
potential between two probe locations 120 and measuring the resultant current, and causing a 
specific current to flow between any two probe locations 120 and measuring the resultant voltage 
drop. 

[0025] Having collected the set of measured resistance values 140, a lumped parameter 
resistance model 160 may be determined for the analysis mesh. The measured values of resistance 
form an NxN matrix [Ry] of measured electrical resistances 140. The overall sheet resistance of the 
conductive film 1 10 may be characterized as a set of lumped parameter resistance models 160. If the 
lumped parameter resistance models 160 are defined between the same probe locations 120 as for 
the measured values of resistance, a matrix of [rij] lumped parameter resistance models 160 is 
defined. A linear transformation, well known in the art, may be developed to define the relationship 
between the lumped parameter resistance models 160 and the measured resistance values 140. This 
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relationship may be defined as a transformation matrix of coefficients [F] based on the physical 
geometries of the conductive sheet and the relative placement of the probe locations 120. Therefore, 
we can write in general that the measured resistance values 140 are a function of the lumped 
parameter resistance models 160. This may be written in matrix form as: 

[Rij] = F[ry] 

i = 1 toN 
j = ltoN 

[0026] As a result there are N equations and N unknowns and one may solve for the 
lumped parameter resistance models 160 using linear algebra and inverting the coefficient matrix to 
arrive at: 

[rij] = F^ [Ry] 

[0027] With a solution for the lumped parameter resistance models 160, we have a 
model to determine resistivity at any point on the conductive film 1 10 encompassed by the probe 
locations 120. This is performed with a weighted average interpolation. 

[0028] The interpolation process occurs by first creating a bounded region of the 
conductive film 1 10 encompassing a selected location 170. The bounded region is defined by 
selecting three probe locations 120 with three measurement lines 150 between the probe locations 
120 defining a triangle encompassing the selected location 170. The lumped parameter resistance 
models 160 for each of the three measurement lines 150 are used in the interpolation. The model 
assumes the resistance per unit length is substantially constant along a measurement line 150 such 
that the resistivity is substantially the same at any point along the measurement line 150. 

[0029] To create the weightings, an orthogonal distance is determined from each of the 
three selected measurement lines 150 to the selected location 170. For example, in FIG. 2, a first 
measurement line 210 between probe locations 120 has a first orthogonal distance 212 (also referred 
to as dza) to the selected location 170. A second measurement line 220 between probe locations 120 
has a second orthogonal distance 222 (also referred to as dzb) to the selected location 170. Finally, a 
third measurement line 230 between probe locations 120 has a third orthogonal distance 232 (also 
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referred to as dzc) to the selected location 170. The weighted contribution of a given measurement 
line 150 is a function of the orthogonal distance from the selected location 170 to the measurement 
line 150 relative to the sum of orthogonal distances for all three measurement lines. For example, 
the weighted contribution for the first measurement line 210 is defined as: 

1-dza /( dza+dzb + dzc) 

[0030] As may be seen and readily appreciated, if the selected location 170 is very close 
to the measurement line 150, the orthogonal distance for that line will be small, resulting in a large 
weighted contribution. 

[0031] The final weighted average is computed as the resistivity (pa, pb, pc) of each 
measurement line (210, 220, 230) multiplied by its corresponding weighted contribution, as defined 
by: 

Pz = pc{ l-dzc /(dza+ dzb + d2c)}+ pb{ 1-dzb /( d2a+ dzb + dzc)} + Pa{ 1-dza /( dza+ dzb + dzc)} 

[0032] A plurality of selected locations 170 may be analyzed to develop a profile of 
resistivity across the entire surface of the conductive film 1 10. Additionally, as described below, the 
profile may be developed for other physical properties that may be correlative to the resistivity 
profile. 

[0033] Selective timing of the measurement and analysis is also contemplated within the 
scope of the invention. It may, for example, be desirable to trigger the measurement and analysis 
after an event which may, for example, be indicated by outputs of one or more sensors. For 
example, if the invention is included on an aircraft wing, it may be desirable to trigger a new 
measurement and develop a new profile after an accelerometer output reaches certain parameters. 
For a storage tank, a new measurement may be taken and a new profile developed when stored 
material within the tank is at a certain level. Pressure and temperature measurements may also be 
used to trigger a measurement and analysis cycle. Alternatively, or in addition to triggering 
responsive to a sensed event, measurements in accordance with the present invention may be taken 
at various times, for example, perhaps periodically, and combined with the profiling to develop a 
time-varying profile of electrical resistivity or other physical property. 
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[0034] For most geometries comprising a relatively large number of probe locations 1 20, 
many different bounding triangles may be defined. As an example, FIG. 3 shows three different 
bounding triangles encompassing the selected location 170. A first triangle is defined by the first 
triangle measurement lines (310, 310', 310") between probe locations PI, P3, and P5; a second 
triangle is defined by the second triangle measurement lines (320, 320', 320") between probe 
locations PI, P2, and P4; and a third triangle is defined by the third triangle measurement lines (330, 
330', 330") between probe locations P2, P3, and P6. Other possible bounding triangles are not 
shown. A more accurate overall weighted average may be possible by combining the weighted 
average from multiple bounding triangles. For example, all possible bounding triangles may be 
combined to arrive at a more accurate overall weighted average. Alternatively, an analysis of most 
likely candidates may be used to identify bounding triangles that may produce the most accurate 
results. For example, an analysis may select only bounding triangles with a combined orthogonal 
distance (i.e., dza+ dzb + dzc) below a predetermined threshold. This would emphasize those 
bounding triangles comprised of measurement lines 150 closest to the selected location 170. 

[0035] An advantage of the present invention is redundancy provided by multiple, 
noncoincident bounding triangles. As explained above, and referring to FIG. 1, with a relatively 
large number of probe locations 120, multiple bounding triangles are likely for any given selected 
location 170. This is an advantage because there may be anomalies in the conductive film 1 10. For 
example, discontinuities may develop due to, as examples only, punctures, tears, cracks, or other 
damage to the surface of the structure. The present invention may be used to locate these 
anomalous areas. However, if a measurement line 150 crosses the anomalous area, the measured 
electrical resistivity value, and as a result the lumped parameter resistance model 160, for that 
measurement line 150 may be inaccurate. The present invention may compensate for possibly 
inaccurate measured resistance values 140 by removing the suspect measurement line 150 from the 
analysis. Typically, other redundant bounding triangles will exist after removal of the suspect 
measurement line 150. These remaining bounding triangles may be used in the weighted average 
analysis to develop a profile for locating the anomalous area. 

[0036] The present invention may be embodied in a variety of physical configurations. 
In the simplest exemplary embodiments, the conductive film 1 10 may be a simple geometric shape, 
such as, for example, the circle of FIG. 1, a square, a triangle, or an ellipse. However, the actual 
two-dimensional shape of the film may take on virtually any shape. For irregular shapes, probe 
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locations 120 may need to be chosen at irregular intervals around the periphery 130 to obtain the 
necessary coverage of the analysis mesh. 

[0037] While the analysis is two dimensional, and the conductive film 1 10 may be 
applied to an essentially two-dimensional structure, the analysis is not limited to two-dimensional 
structures. Rather, when used on a three-dimensional structure the analysis is of the surface 
properties of the structure, as opposed to the volumetric properties of the structure. The film may be 
applied to curved surfaces or across a plurality of surfaces comprising the three-dimensional 
structure. Thus, the conductive film 1 10 may be attached to various structures such as, by way of 
example only, storage vessels, ship hulls, aircraft wings, spacecraft wings, turbine blades, body 
armor on military combat tanks, personal body armor, and vehicle axles. FIGS. 4A-4C illustrate one 
exemplary three-dimensional structure 410. FIG. 4A illustrates a rectangular conductive film 1 10 
including probe locations 120 and measurement lines 150. FIG. 4B illustrates a cylindrical storage 
vessel 410. FIG. 4C illustrates the conductive film 1 10 attached around the perimeter of the 
cylindrical storage vessel 410. As can be seen, the probe locations 120 are still around the periphery 
130 of the cylindrical storage vessel 410 at points along the top circular surface and the bottom 
circular surface. This configuration allows detection of anomalies around the cylindrical portion of 
the cylindrical storage vessel 410. Additionally, depending on the type of conductive film 1 10 used, 
this configuration may enable detection and profiling of other physical properties such as stress 
characteristics, thermal characteristics, and photosensitivity. 

[0038] It should be noted that the analysis may not be as effective when using periphery 
measurements on a structure or film with a very large aspect ratio. For example, the cylindrical 
storage vessel illustrated in FIGS. 4B and 4C, has a large aspect ration when the cylindrical storage 
vessel is very tall, with a small circumference. If probe locations are only at the top and bottom 
circles, the analysis triangles become significantly extended in one direction and small along at least 
one of the triangle legs creating diminished resolution along the long axis. If intermediate probe 
locations may be placed along the long axis of the structure or film, this large aspect ratio disparity 
may be alleviated. 

[0039] A variety of conductive films 1 10 may be employed in the present invention. As 
long as the conductive film 1 10 has adequate electrical conductivity for the analysis described 
above, the conductive film 110 may be adapted to be sensitive to other physical properties. For 
example, the conductive film 110 may be made of a metal having a correlation between resistivity 
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and deformation, similar to that of metals used in strain gauges. A conductive film 1 10 having this 
correlation to deformation may be used to profile stresses across the conductive film 1 10 and 
surface of the structure attached thereto. As another example, the conductive film 1 10 may comprise 
a metal film similar to that used for metal film temperiature transducers, such as platinum, enabling 
profiling of temperature across the conductive film 1 10 and surface of the structure attached thereto. 
A temperature profile may also be extrapolated to related parameters such as thermal transmissivity 
of portions of the underlying structure Yet another example is a photosensitive material, such as for 
example, cadmium sulfide. Using a photosensitive material enables profiling various intensities of 
light impinging on the conductive film 1 10. Additionally, a photosensitive material may be targeted 
at specific radiation wavelengths. 

[0040] The mode of attaching the conductive film 1 10 to the surface of the structure also 
may vary. Attachment methods may vary depending on the material used for the conductive film 
1 10 and physical properties of interest. For example, if a deformation of stress at various points on 
a structure is desired, the conductive film 1 10 may be applied using an adhesive that enables the 
conductive film 110 to contiguously deform, compress, or stretch with the underlying surface of the 
structure. On the other hand, if surface temperature of the structure is the desired property, the 
conductive film 1 10 may be attached in a different manner. For thermal applications, it is more 
important to ensure adequate thermal conductivity between the material to be measured and the 
conductive film 1 10, rather than the distributed physical attachment required for stress 
measurements. 

[0041] The conductive film 110 may not necessarily comprise a prefabricated sheet 
adhered to the structure in some fashion. Instead, the conductive film 1 10 may be applied to the 
structure using a method such as spraying on the film, the spraying technique varying with the 
material of the film. For example, a metallic coating may be thermally sprayed on a surface. Thus, a 
powder of a material such as alumina or copper may be sprayed at a high temperature, using so- 
called plasma spraying techniques, onto a structure creating a thin, continuous, and evenly 
distributed conductive film 1 10. For extremely fine control of film properties, chemical vapor 
deposition (CVD), plasma enhanced chemical vapor deposition (PECVD), or even atomic layer 
deposition (ALD) may be employed to deposit a variety of materials. 

[0042] As shown in FIG. 5, the present invention also includes a system configured for 
determining surface properties of a structure bearing the conductive film 1 10 over the structure's 
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surface. This system includes a plurality of probes 510 adapted for measuring electrical resistance 
140. Each of the plurality of probes 510 connect to a signal controller 520. The signal controller 520 
typically may be an analog multiplexer configured for selecting any two of the plurality of probes 
510 for making a resistance measurement across a selected measurement line 150 (see, for example, 
FIG. 1). After making a resistance measurement, the analog multiplexer may be switched to a 
different pair of probes 510 to make a measurement on a different measurement line 150. A signal 
sampler 530, connected to the signal controller 520, may be used to sample a value on the currently 
selected pair of probes 510 and convert the sample from an analog signal to a digital signal. A 
suitably programmed processor 540, connected to the signal sampler 530, may be used to receive 
the digital signal and performs the analysis described above. The processor may be any computer, 
microcontroller, microprocessor, digital signal processor, or custom circuit, configured for 
performing the required analysis. 

[0043] Although this invention has been described with reference to particular 
embodiments, the invention is not limited to these described embodiments. Rather, the invention is 
limited only by the appended claims, which include within their scope all equivalent devices or 
methods that operate according to the principles of the invention as described. 
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